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Abstract. We write equations of motion for density variables that are equivalent 
to Newtons equations. We then propose a set of trial equations parameterised by 
two unknown functions to describe the exact equations. These are chosen to best fit 
the exact Newtonian equations. Following established ideas, we choose to separate 
these trial functions into a set representing integrable motions of density waves, and 
a set containing all effects of non-integrability. The density waves are found to have 
the dispersion of sound waves, and this ensures that the interactions between the 
independent waves are minimised. Furthermore, it transpires that the static structure 
factor is fixed by this minimum condition to be the solution of the Yvon-Born- 
Green (YBG) equation. The residual interactions between density waves are explicitly 
isolated in their Newtonian representation and expanded by choosing the dominant 
objects in the phase space of the system, that can be represented by a dissipative term 
with memory and a random noise. This provides a mapping between deterministic and 
stochastic dynamics. Imposing the Fluctuation-Dissipation Theorem (FDT) allows us 
to calculate the memory kernel. We write exactly the expression for it, following two 
different routes, i.e. using explicitly Newtons equations, or instead, their implicit form, 
that must be projected onto density pairs, as in the development of the well-established 
Mode Coupling Theory (MCT). We compare these two ways of proceeding, showing the 
necessity to enforce a new equation of constraint for the two schemes to be consistent. 
Thus, while in the first 'Newtonian' representation a simple gaussian approximation 
for the random process leads easily to the Mean Spherical Approximation (MSA) for 
the statics and to MCT for the dynamics of the system, in the second case higher levels 
of approximation are required to have a fully consistent theory. 



1. Introduction 

The field of equilibrium statistical mechanics was founded on the simple but remarkable 
proposition that motions generated by Newtonian or Hamiltonian mechanics have very 
simple statistical properties if there are sufficiently large numbers of interacting particles. 
Specifically, the earliest developments of statistical mechanics presupposed the so-called 
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theorem of molecular chaos which assumed that all parts of phase-space are equally 
probable [|TJ. This simple idea is at the basis of the ergodic theorems, which state the 
equality of time and ensemble averages, the latter having an explicit representation in 
the Boltzmann probability distribution. Subsequent developments further strengthened 
the foundations of statistical mechanics by weakening the molecular chaos requirement, 
and the hypothesis of equal probability for all parts of phase space, while retaining 
the formalism of statistical mechanics. Other related developments allow to study the 
dynamics of many particle systems using some sort of stochastic process that mimics 
aspects of the chaotic motions of Newtonian mechanics. Indeed, various advances have 
provided a firmer basis for relating Newtonian mechanical averages to the statistical 
properties of stochastic averages |2|, |3| . 

Somewhat deeper problems arise when one wants to study those systems where the 
parts of phase space explored over time are limited, and limited in what is perhaps no 
longer a simple manner. Already approaching an equilibrium phase transition, a system 
becomes strongly non-ergodic. Indeed, when a phase transition occurs, by the breaking 
of a symmetry for example, it is no longer appropriate to average over all parts of the 
phase space in the probability averages of Boltzmann. Instead we average only over 
those parts close to the lowest, most probable energy state J4j. In such simple cases the 
Boltzmann machinery, applied to a limited portion of phase space, can still be applied. 
However, when the ergodicity is lost in the system in a less well defined manner it may 
not be possible to find any straightforward generalisation of the Boltzmann distribution, 
and one is forced to reconsider the possibility to perform time averages with respect to 
another type of averaging. The provision of a well defined mechanism to do this is, in 
the arena of spin glasses, one of the substantial achievements of that field ||. 

Where dealing with particles, whether they be atoms, molecules, or colloidal 
particles, we know at least that Newtonian dynamics, or some equivalent stochastic 
process, remains faithful to the physical picture, and averaging with respect to such a 
dynamics would be satisfactory. Examples of systems where these questions arise include 
structural glasses ||, the glass transition f7j, and many other cases where amorphous 
materials such as gels [§, || |10[] , composites etc. are formed, there being no prominent 
state to dominate the average, and render the Boltzmann procedure acceptable. In such 
cases one of the few approaches able to describe the dynamics of the system is the ideal 
Mode Coupling Theory (MCT) |]lTJ, |12|]. In fact, as we shall see, one way of viewing 



MCT is as an approximate method of performing Newtonian averages [13]. In this sense 



MCT is an important exemplar of a new way of thinking about calculating statistical 
dynamical averages, without depending on any form of ergodic hypothesis, or on any 
assumption that the long time probability distribution is the Boltzmann one. 

In this paper, we shall present a formalism that connects Newtonian averages taken 
with respect to initial positions and momenta to stochastic averages taken with respect 
to noise distributions. This formalism has much in common with the Zwanzig-Mori 
formalism 0, but we here emphasise as much as possible the connection and fidelity 
to Newtonian averages, and the possibility to explicitly realise these averages within 
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various approximations. 

An important aspect of the approach lies in its transparency. This allows the 
possibility of an easy generalisation, for example, to consider other types of collective 
variables, other than simply density waves, as well as to represent non-integrable parts 
of motion of the system not only by linear dissipation and random noise, as we will 
discuss here. We will also attempt, where possible, to probe the physical meaning 
of the various steps, and above all the limitations of the approximations made. One 
interesting consequence of our approach is that ideal MCT emerges as a consistent and 



natural first 'mean-field' type approximation [|13| |14| , |15|1 to our theoretical framework, 
thereby providing an alternative route to the MCT as a theory of dynamics, consistent 
with the Random Phase Approximation (RPA) of statics. In this regard, we can now 
consider MCT to be exact in the weak coupling limit. 

Thus, we begin by writing equations of motion for the density variables that are 
equivalent to Newtons equations. We then propose a set of trial equations, parameterised 
by a set of unknown functions, to describe these exact equations [[U|. These parameter 
functions will have to be chosen to best fit the exact Newtonian equations. We choose to 
separate the set of trial functions in the trial equations of motion into a set representing 
integrable motions of density waves (the dispersions), and a set containing all the 
remaining effects of non-integrability. In so doing the density waves are seen to have 
the dispersion of sound waves, and this choice ensures that the interactions between the 
independent density waves is minimised. Furthermore, it transpires that the structure 
factor, and thereby the direct correlation function, is fixed by this minimum condition 
to be the solution of the Yvon-Born- Green (YBG) equation 0, [T7|. As part of this 
process, residual interactions between the density waves are then explicitly isolated in 
their Newtonian representation. We therefore have an explicit Newtonian representation 
of the interaction force, and find a way of representing this force by expanding the trial 
function. 

Evidently it is possible to expand these interactions in a set of functions that are best 
chosen to represent the dominant objects in the phase space of the system that provide 
for loss of the pure oscillatory motions of density waves. The key question therefore is the 
nature of the underlying phase space in the density representation for that state of the 
system we seek to describe. It is emphasised that these interaction parameters represent 
all the non-integrable contributions of density waves, rather than particles. Thus, the 
phase relations between density waves that are implied by the existence of particles, as 
well as the interactions of these particles and the resulting dissipation and other type 
of non-integrability of the underlying Newtonian system must be accommodated by the 
trial functions. To make these ideas more concrete we will illustrate the problem with 
a simple example of how the non-integrable parts of phase space can be represented 
by a dissipation term with memory, as well as a correlated noise to represent the more 
rapid fluctuations of the density caused by effects like chaotic motions. This leads to a 
generalised Langevin equation where the distribution of noise and the memory kernel 
of the dissipative term are to be determined to reflect the statistical properties of phase 
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space. 

It should be emphasised that, in choosing to approximate phase space statistics 
partly with deterministic and partly with stochastic terms, we will always naturally 
arrive at a stochastic process for the trial density function equation. Therefore the idea 
of 'best' representing Newtonian phases space motion can only be implemented in some 
average sense, rather than point-wise in particle position and momentum. This will 
mean that we must relate, and equate, general averages over the stochastic variables to 
those over the Newtonian space. The equality of these two averages will, at each order, 
define the 'goodness' of the approximation. These issues are discussed in section [5|. 

Returning to the simple case where the trial forces or interactions are described by 
a noise and a memory kernel, to produce a conventional generalised Langevin equation 
it will be sufficient to apply the Fluctuation Dissipation Theorem of type two (FDT) 
T8fl to ensure that the system can reach equilibrium in the long time limit. It is natural 



to ask whether the use only of these two terms in the trial forces is sufficient to represent 
Newtonian dynamics, and we shall discuss the question at some length in section |5.2| . 

We have earlier alluded to the matter of ensuring that averages over density 
correlators are the same in Newtonian and stochastic spaces. Since this involves a 
hierarchical process, there is the practical need to truncate the program at some point. 
One practical possibility, outlined in section |573| , is to ensure equality of the relevant 
correlators in the two spaces, and in particular to ensure that the stochastic part of the 
force is correlated properly with the density variables, again to some appropriate order. 
These considerations then lead us to a hierarchy of equilibrium correlation function 
equations, the first level being the YBG equation discussed above. 

Following this, in section we turn to the evaluation of the exact memory kernel, 
using the two alternative representations of the stochastic part of the force. This 
culminates in two exact equations ([45D and (|46|) , that represent the central results of 
this paper. While these equations are exact, they are not closed. Therefore, in sections 
|7|-|] we discuss various aspects of closure strategies. The most elegant closure found 
to date involves assuming a Gaussian distribution for the noise, using the Newtonian 
represntation, and this leads us (section [5]) directly to the known ideal MCT for 
supercooled liquids. 

It is important to note that, while our formulation and presentation are different 
from earlier treatments, in many aspects there are parallels and links to previous work. 
Thus, in section |]we show how to deal with the implicit representation of the force, and 
point out that only pairs of density variables are required for its exact specification. We 
also, in section |8.2| , note that the Gaussian approximation can only be consistent up to 
first order in density variables, and also that the stochastic force cannot be consistently 
represented at the level of density pairs if it is gaussian distributed. Next, in section [H7^, 
we begin the process of seeking consistent approximations that permit the representation 
of the stochastic force at the level of pairs of density variables. As an example, we 
illustrate the Singwi-Tosi-Land-Sjolander (STLS) [|T^] closure is introduced as a step in 
this direction, but it also fails to be completely satisfactory. 
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Finally in section |9], by extracting the main elements of the traditional derivation 
of MCT via projection operators, we will attempt to clarify the assumptions made in 
our derivation of MCT and in the traditional treatment [O . 

2. Density Equations of Motion 

We begin with the definition of Newton's equations, 

1 9V 

where Tj(t) are the positions of particles j — 1, N and m is the particle mass, and we 
suppose that the total inter-particle potential V(r 1; r n ) can be expressed as a sum of 
pair potentials, 

V = \Y, v{\v rn {t)-v m ,{t)\) (2) 

m,m' 

The density of the N particles located at positions ij is p(r, t) = J2f=i <K r — r j(^)) an d 
its Fourier transform, which we name the density variable, is, 

N 

Pk(t)=E^ W (3) 
and the structure factor of the system is defined as the density-density correlation, 

= l(p_ k (0)p k (t)> (4) 

which, for equal times, corresponds to the static structure factor Sk = (|pk| 2 )- 

We, thus, write the equations of motion corresponding to these variables using the 

second time derivative of (|3p, 

p k (t) = £i(k • f s {t))^'^ - £(k • r,(t))V k ^ (5) 
j j 
So far these equations are only formal equalities, but we can now eliminate the second 
time derivative of the positions using Newton's equations ([!]). To do this, we Fourier 
analyse equation (|j) to obtain, 

V = ^EE vve-^V-^'W = ^E vA&PK - N5 k ,, Q ) (6) 

m,m' k' k' 

where V is the total volume of the system and we have assumed that the system is 
isotropic, and thus, the potential depends only on the modulus of the wave-vector. 
The derivative of the potential in Fourier space ([]) can be written as, 

so that the equations of motion for the density variables (0) now become, 

p k (t) = - y(k ■ rJtW^j® - — YYv k ,(k • k ') e -*'-(r,(*)-r m (t)) e ik.r 3 W 

= - E(k • ijit))****'® - E^(k ■ k)p k -k'(t)Pk'(t). (8) 

■j m v k' 
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Note that equations (||) constitute an exact set of equations, that are equivalent to 
Newton's equations, but written in terms of the density variables. Clearly the first 
term represents the kinetic part of the motion, which is not simply expressable in terms 
of density variables, while the second term, originating from the interaction potential, 
indicates a coupling of the density pairs. 

We note here that, in principle, to describe the underlying Newtonian dynamics 
correctly, i.e. to include sufficient information in the theory that would mimic positions 
and velocities of the N particles, one would need to consider many fields. In particular 
the most general ones are usually considered the hydrodynamic variables of longitudinal 
and transverse currents, along with a local temperature or entropy variable [p0| . Until 
now, mode coupling type theories have been based on longitudinal current, and therefore 
only the density degree of freedom, and this is the approach we will also follow. 
Extensions involving other fields are more complex, but, due to the transparency of 
the approach, it should be quite straightforward at least to write the correct equations 
of motion, as done in this and the next sections for the density variables. 



3. Trial equations and Residual Force between Elementary Excitations 

Equations @, being equivalent to Newton's equations, are intractable as presently 
written. They must therefore be approximated in some manner. To do this we consider 
the possibility to isolate any simple degrees of freedom that can be easily treated, 
especially if the real system can be considered as some sort of perturbation around 
them. Natural candidates for this treatment are all integrable motions in the phase 
space. We will consider as fundamental collective motions for describing the system and 
the essential mechanisms by which its fluctuations relax the density waves. 

Thus, we propose to write equations (j§) in a simple 'trial' form Indeed, we 
consider the equations of motion, for density waves that have some interaction force 
between them, to be of the form, 

Pk(t) + fi k p k (t) = A(t) (9) 

Here, we have isolated the linear term in Pk(t), that represents the elementary excitations 
of the system, and we have indicated by JF k (t) the interaction force between them, this 
being a functional of the frequencies Cl^ that have to be determined, as well as any other 
parameters to represent the forces. These elementary excitations can be interpreted as 
a sort of 'phonons'. Of course, these would be oscillating perfectly with the frequencies 
fik if there was no interaction force between them. Thus, they represent integrable 
motions and we also assume that, for the purpose of describing a liquid, these are the 
only collective motions that need to be taken into account. The residual interaction 
force between these 'phonons', Fk(t), is the responsible for the deviation of the density 
waves from a perfect oscillatory behaviour. The smaller it is, the simpler the solutions 
Pk(t) of equations (0) will be. Generally, in a liquid, this interaction is not negligible, 
as it would be in a crystal. 
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Now, we apply a variational principle to minimise this interaction [16|, |21|. Our 
variational parameters are the dispersions f2 k . We require, 



<9<l^k| 2 > 



(10) 

and thereby calculate the optimal dispersion f2 k , that satisfies this condition and thus 
results in the least residual interaction between density waves. 

Here the averages are performed over initial conditions of positions and velocities, 
{r n (0), r n (0)}, or equivalently of density variables and their velocities, {p k (0), pk(0)}, 
via the definition (H). Therefore we have indicated them with the notation (• • to 
denote that they are Newtonian averages. We have, 

d(\A\ 2 )« d(|Pk(t) + lWt)| 2 )^ 

<9f2 k <9f2 k 
so that the solution is, 



(P-k(Pk + ^kPk))jV+((p'-k + ^kP-k)Pk)jV = (11) 



^k = jj^ (12) 



This is a general result. Indeed, it was shown by Zwanzig jn| to be valid for the 
frequencies associated with any elementary excitations in a fluid, in our case p k (t). By 
using the definitions of density variables (|3|), one easily obtains, 

fik = /V (13) 



The same result (13) can also be easily found by evaluating the short-time limit of 
equations (H), after multiplying them for p_ k (0) and performing the thermodynamical 
average |3|, |l2fl . The particular derivation we have presented here is instructive in 
showing how the collective motions are being removed from further consideration in 
the dynamics of the system, as all of the attention from now on can be focused on the 
residual interactions between them. We can now calculate this explicitly. By adding to 
both sides of equations (0) the term f4p k (£), we have an explicit form for the minimum 
residual force JF k (t) = minjF k (t). Thus, 

^k(t) = ^Pk(t) - E( k ■ *i(t)) 2 e*- r ' (t) - 4? EM k ' k')Pk-k'(t)/v(t). (14) 

3 k' 

Clearly at some point we will have to represent this force in a simpler manner, if we are 
to make progress. 



4. The Yvon-Born-Green Equation 

It is of interest to note the fact that the best excitations (in the sense of minimising 
residual interactions between them) also implies some conditions on the interaction 
force. Indeed, if we refer back to the minimisation of equation (0), we can see that 
there an orthogonality condition is implied, i.e. 

(p-k(*)^k(*)>* = (15) 
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This is the statement that the instantaneous force at every time must be orthogonal to 
the density waves themselves. This is an exact condition that should be enforced in any 
type of 'trial' equations of motion developed to describe Newton's equations. Thus, it 
must also be applied to our equations, with minimised residual force, 

/3k (t) + ^kPk(t) = ^k(t) (16) 

Of course, this equation can be used to determine the equilibrium correlation functions. 
Indeed, using the explicit expression of the residual force ([14]), the orthogonality 
condition may be written out explicitly as, 

fi k (|Pk|%-(E( k ' r,W)V kr ^-^ 5> k ,(k • k')(p- k (0Pk-k'(*)Pk'(0)^ = 0(17) 

j k' 

Now we note that the averages involved here are equilibrium averages, and as such they 
involve nothing more than the Boltzmann distribution. Thus, the averages of positions 
and momenta (or velocities) factor and can be carried out separately to yield, 

k 2 / 1 \ 1 



, , ^)NS k = — -E^(k-k')(p- k (t)Pk^WPk'(t))^ (18) 
pm \ ok/ rnv k , 

and defining the direct correlation function, Ck = \ (l — where n = N/V is the 
number density of the system, we find, 

c fc = -I3v k - ^ J2 u k'(k • k')(p- k (*)Pk-k'(*)Pk'(*))Af (19) 

ly ^ k'^k 

where we have extracted from the sum the term corresponding to k = k' since it 
corresponds to simply —(3v k - Condition flT9"| ) constitutes a set of equations for the 



equilibrium fluid. It can also be shown W% that Eq.|19] is simply the representation 



in Fourier space of the YBG equation. Thus, equation ( |TP| ) gives a relation for the 
pair (static structure factor) and triplet density correlation functions in terms of the 
potential of the system. For example, if one could determine the static structure factor 
by ensuring this condition, one would have also determined the set of frequencies Q k 
that appear in the trial equations (p~6|) . 

5. Transition from a Deterministic to a Stochastic Dynamics 

As we commented under equation ([14]), it is necessary to approximate in some way the 
residual force to make further progress. Broadly speaking, it has always been assumed 
that such forces can be separated into a slowly varying and a rapidly varying part 



|2|, |3], ^2j. If we decide to represent the rapidly varying part as some type of noise, then 
the trial equations of motion for the density variables will become stochastic equations. 

In this section we will examine the transition from a deterministic to a stochastic 
dynamics in its various aspects. Firstly, in subsection |5.1] , we will propose a simple form 
for the residual minimised interactions T^t) between density waves, basing ourselves on 
standard mathematical formalisms ||. This will allow us to write the density equations 
of motion (||) in the form of a set of generalised Langevin equations. Then, in subsection 
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572 we will focus on the mathematical constraints that must be imposed to our stochastic 



process to faithfully reproduce its 'true' Newtonian dynamics. Finally, in subsection ^3 
we propose a practical way of implementing constraints between the Newtonian and the 
stochastic space. 

5.1. The Physical Picture 

So far, we have written a set of trial equations (|T6|) , where all of the collective parts 
of the motion have been isolated and, at the same time, the interactions or effective 
forces between these collective motions have been minimised. The outcome is a set of 
collective motion frequencies fllTf ) that are related to the structure factor, and that this 
structure factor is the solution of the YBG equation. 



Now, the effective interactions, or trial forces, in equation (0), must be approximated 
to accommodate the residual physical interactions, but preserving the mathematical 
constraints, such as ([15]) discussed above, between them. A central observation is that, 
in the absence of the coupling terms T^t), the collective motions would have infinite 
lifetimes. On the contrary, in a real system, depending on its details, we have different 
types of effects, tending either to cause the collective phonon-type motions to decay or 
to make them more long lived. Indeed, the same system, for different temperatures 
or densities, may exhibit such different tendencies, so that it will be important to 
design trial forces that are flexible enough to accommodate the particular phenomena 
of interest, as for example, the glass transition. 

Let us first deal with the issue of particle versus wave motion. Of course, if we 
are able to solve Newton's equations exactly in either their familiar particle or their 
density representation there will be no difference. That is, we may represent all of the 
phenomena equally well in terms of particles, or density waves, providing every step is 
exact. Inevitably approximations will be introduced and then the representation that 
we choose will be important. For smaller wave- vectors we know that a substantial 
part of the motion is collective. For larger wave-vectors we expect that aspects of the 
independent and uncoordinated nature of the single-particle motion will be dominant. 
It should be noted that at long wavelengths we are dealing with waves where many 
particle collisions occur and where the longer density waves are locally in equilibrium. 
This is also reflected in the fact that we seek to model the collective oscillations using 
frequencies that involve the equilibrium structure factor. Therefore it is only part of the 
collisional spectrum that contributes to the decay of the density variables, the remainder 
having contributed to local equilibration of the density waves. 

We now turn to the question of how the remaining effects of finite lifetimes of these 
waves can be described. Also, we recall that it will be important to retain at least 
some information on the relative phases of the density waves, so as to have an adequate 
description of particle, rather than wave, properties. 

We assume, as usual, that the finite lifetime arises from correlated phenomena, as 
well as more chaotic effects. It is thus conventional to separate the residual interactions 
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into two types of contribution, following the ideas of the Zwanzig-Mori formalism 0. 
The first term involves a first derivative in density, so it will clearly cause a decay of 
the oscillating density waves. However, by using a memory kernel, we can parameterise 
some aspects of the history of the system. The second term is a general random noise 
that originates a more 'chaotic' behaviour for the system. Of course, such a separation 
is not a final and complete one, but it is conventional, and we will accept it for the 
current purposes. 

By this statement we can represent our residual force (fU ) in the most general form 
using these contributions to non-integrability, 

^k(t) = ~ f 7k(t - t')p\{t')dt' + 77 k (t) (20) 
Jo 

where 7k (t) is the memory function of the system and rj^t) is the noise. The form (ED) 



for the residual force must, of course, satisfy the orthogonality condition (|T5|). 

We now return briefly to the effects of having a memory kernel in equation (|20|). 
Clearly, by selecting various forms for 7k (t) we are able to control the means by which 
dissipation is expressed in the system. Thus, short-time effects that would, for a normal 
liquid, lead to more rapid decay of the waves, can be incorporated. Also, more extensive 
collective effects involving many particles can be incorporated, in an averaged manner. 
Finally, again in some averaged manner, we can incorporate those aspects of phase 
coherences of different density waves, that are necessary to retain a minimal picture of 
particle motions. Indeed, potentially we are using the memory kernel to accommodate 
a number of different physical effects. 

By this expansion, we have made a transition from a deterministic description of the 
system, based on equations ( |16|) and ([14]) , to a stochastic one, where the full equations 



of motions are expressed through the generalised Langevin equations, 



Pk(t) + ^fcPk(t) + / 7k(* - t')p k {t')dt' = r) k (t) (21) 



o 



with f2k given by (|T^)and 7k (t) still to be determined. Introducing the normalised 
density correlator, defined as $&(£) = Sk(t)/Sk, the exact equations of motion clearly 



$k(t) + fik$k(t) + / 7k(t - t')K{t')dt' = (22) 
Jo 

where the causality relation (p„k(0)r?k(£)) = for t > must hold. 

In our system the fluctuations, represented by r/k{t), and the dissipations, expressed 
by 7k(£), are n °t independent, and, in particular, they must satisfy the fluctuation- 
dissipation theorem (FDT), or else no Boltzmann equilibrium distribution would be 
reached by the Langevin equation implied by equations (|2~lD . Thus, we have, 

(Ipkl 2 ) =7k(t -° (23) 

This relation fixes the memory function. It can be easily shown that this relation is 
directly implied by our generalised Langevin equations fl2"T| ) , with the only assumption of 
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equilibrium. To describe out-of-equilibrium dynamics, a generalization of this formalism, 
and consequently of FDT, is necessary, as we discussed in p3| . 

We now want to clarify the question of time-translational invariance of our 
stochastic process ([H|). Indeed, there is an explicit dependence on the arbitrary initial 
time t — 0, at which memory effects start to be taken into account, due to the crucial 
choice for the residual force (|2"0|). This is done for obvious physical reasons. If we 
consider another Langevin process of the form, 

p k + n kPu (t) + [ 7k (t - t') P \(t')dt' = f) k {t) (24) 

J — oo 

which is now clearly time-translational invariant, we have the following relation between 



our true random noise r/^t) of equations (^j]), and this new, auxiliary, noise r/^t) H 

Vk{t) = fj k {t) - f° dt' lk {t - t')p k {t') (25) 

J — oo 

it is sufficient to ask our non-time-translational invariant noise rj k (t) to obey the initial 
conditions, 

(p-k(0)77 k (t)> = t > 

(p-k(OMt)) = t > (26) 

as well as that its autocorrelation function is the same as the one for ffc(t). Of course, 
we require the memory kernel, 

(77-k(0)»fc(t)) = (r?-k(*i)7?k(ti + t)) = 7kW (27) 

to ensure that the two stochastic processes (|21| ) and (|24"D are completely equivalent |T8|] . 
Thus, causality is also correctly built in our equations (^Tj). 

This means that, though the random noise rj^t) itself is not time-translational 
invariant, the whole process, described by equations (0), is. In particular, we have that 
the residual force J-^if) (p0|), which is the sum of the random noise plus the conventional 
memory term, starting from t = 0, does not depend on the initial time. Indeed, it is 
important to clarify this aspect, as the orthogonality condition (|1|) between the density 
and the force T at equal times, implies, 

(p-k(t)^k(t)) = (p-k(O)^k(O)) = (p-k(OMO)) = (28) 

since at the initial time, we have J^k(0) = ^k(O). But we cannot simply say that also 
(p_k(t)?7k(^)) is zero, because of the non- validity of the time-translational invariance for 
the noise ffk(t). In fact, for any time t, the orthogonality condition just implies, 

(p-k(*K(0> = /V(* - t')(p- k (t)p k (t')) (29) 
Jo 

and, in particular, each of these terms is zero in the trivial case of a ^-function memory 
kernel, for which, of course, the equations motion become local in time. 
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5.2. The Mathematical Constraints 

In this section we will discuss the mathematical definition of our stochastic process, and 
the implied constraints. This will shed some light into how good our Langevin process 
([2~T1) is in reproducing the Newtonian process in equations (||). Thus, the stochastic 
process r]k(t) can be defined via the generating function, 

{tj) ) {Vk} = (e*^ fe)+ /oS fe-*')^ (i ' )d '' ] ) {rn(0)irn(0)} (30) 

where the average on the lhs corresponds now to an average over the noise probability 
distribution. Equation ([3C]) is, in principle, a full definition of the random process, 
because a stochastic process is completely defined by all its moments. Indeed, every 
correlation function of the noise, considered as some moment of the distribution, could 
be calculated by deriving it with respect to the auxiliary parameters and it is thereby 
related directly to an average with respect to Newton's equations. These relations are 
only just formal, but they do place a clear set of constraints on the stochastic process 
that we choose to mimic Netwonian dynamics. 

Besides this relation between the stochastic process and Newton's equations, there 
is also, via equation fl2~T|), a mapping between the density variables of the Langevin 
process and the noise. Indeed, Laplace transforming equations fl2T|) and solving with 
respect to the density, we have, 

~ , \_ H* + 7k(z)]Pk(0)+Pk(0)+i/k(z) ,.,-n 
-z l - izj k (z) + ll k 

where Pk(0) and Pk(0) are the initial conditions in time of our variables, and we have 
used the definition of Laplace transform, 

POO 

p k {z) = / dte izt p k {t) (32) 
Jo 

z being a complex number. Thus, now averages over density variables in the 
deterministic language can be replaced by averages over the noise r/^it) in the stochastic 
one, or their equivalent density averages. 

Equation (|23|) defines the autocorrelation function of the noise distribution, but to give 
a complete definition of the stochastic process we also need to define the mean value 
of the distribution. To do so, we can use the mapping with density variables, which 
implies the further condition (^(t)) = 0. 

Now, in principle, to create the noise distribution we should use the definition of 
the stochastic process, via the generating function in ([3(]), and by taking derivatives of 
this, calculate the noise distribution moments satisfying this equality up to any desired 
accuracy. From these constraints, then, one can determine the density correlators. 
The level at which one truncates this program, via some closure, will be the level of 
approximation of the theory. Evidently, if one followed this program, all averages over 
any function of density variables taken in the Newtonian space could be explicitly written 
in terms of averages over density variables in the stochastic averages. 

However, in practice we are not interested in the correlators of arbitrarily high 
order, nor are we interested in arbitrary functions of density variables. Indeed, what we 
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still need to determine to have a complete theory is the memory kernel 7k (t), which must 
satisfy the FDT relation (^). We will illustrate this in the following of the chapter. 
Nonetheless, for the present purpose, it is interesting to anticipate that the result for the 
memory will depend only on 2-, 3-, and 4-point density correlators for unequal times, 
and thus these are the only averages that need to be equal in the two spaces, Newtonian 
and stochastic respectively. Therefore, we should focus most attention only on these 
correlators. 

We can write explicitly the implied equalities for 2- and 4-point correlation functions 
at unequal times for the noise, derived from (|30|) , i.e. 

rt 



(^WWOW = (^(0) \F* 2 (t) + J q dt'^t-t')^' 



{r n (0),r„(0)} 



(33) 



(77_ kl (0) 77_ k2 (0) 77 k3 ) ^ (* ) > {rMc} = (V_ kl (0)^_ k2 (0) \j^(t) + J*dt' lk3 (t-t')p ks (t' 

?U (t) + f dt' lu (t - t')p k4 (t')l ) (34) 

L JO J / {r n (0),r n (0)} 

These equations define the stochastic averages of the correlations, or what is equivalent, 
the definition of the stochastic process in terms of the Newtonian dynamics, to the 
appropriate order. But we can further note that, using the definition of ^(t) in terms 
of the density variables ( [14] ) and substituting into this the expression for the densities 
in terms of noise (|3~i"l), we obtain a hierarchy of constraints for the noise distributions. 
Indeed, without writing explicitly all of these contraints, it is easy to see that equation 
(j33[) implies a relationship for (7777) in terms of (777777) and (77777777). Similarly, equation 
(j34f) implies relations between (77777777) and higher order correlations of the noise up to 
order (7777777777777777) . 

These constraints should be enforced for any implementation of Mori-type formalisms, 
with respect to the appropriate Newtonian dynamics. Then, practically, all closures of 
the theory involve violation of these constraints at some order of noise moments. 



5.3. Equivalent Representations of the Random Force and Practical Implementation of 
Constraints 

We have seen so far that the right constraints to implement for the stochastic process 
to represent faithfully the Newtonian process follow directly from the definition of the 
random noise fl30|). More practically, we have written explicitly the type of constraints 
that will be requested for our simple theory, based on the expansion (p0|) for the residual 
interactions T^t), in (^) and fl3"4|). Now, if one correctly implemented the constraints 
of type fl33"D and (|34]), as well the similar ones for the triplet correlation averages that 
will appear in the expression for the memory kernel in the next section, the transition 
from Newtonian to stochastic averages would be successfully achieved for our formalism. 
Thus, in particular, one could consider the two sides of the exact density equations (|5|) 
to be equal when averaged with respect to any order of the density variables using the 
Newtonian phase space, or via the relation (|3~1~D , to any order of the noise variables. 
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However, the execution of this program, i.e. the exact implementation of (^) and 
(SI]), is generally difficult to pursue, and it may be of interest to ensure less extensive 
conditions, but ones which address the most important physical effects. The fact that the 
two sides of equations (^) should be equal with respect to the averaging process, to some 
chosen order of the density variables, at least ensures that the underlying Newtonian 
dynamics is reasonably reflected in the stochastic process. Thus, considering equation 
( pOl) , it is clear that the random force rj^t) can be explicitly written as, 

rtk(t) = Mt) + f 7k(t - t')p k (t')dt' (35) 
Jo 



but the residual force can be expressed either in its explicit form, given in fll4D, or in its 
implicit one, given by flTo]). Thus, we have the two equivalent representations for the 
noise, 

^ /) (t)=fi fc p k (t)-^(k-r i (t)) 2 e lk -^W ^^ k - k ')/^(*)/^(*) + A* , 7k(t-0Pk(* , X36) 

and, 

V ( "\t) = W) + ^Pk(t) + f dt' lk (t - t')p k (t') (37) 

Jo 

For convenience of notation, we will refer from now on to eq. (0) as representation I of 
the random force (or explicit representation), and to eq. ([37]) as representation II ( or 
implicit representation). The requirement that the two are equal, which is to say that 
the decomposition of the residual force (|20|) is valid, is clearly equivalent to consider the 
equality between the two sides of equations (§), since the integral term, that we define 
for simplicity, 

A k (t) = f\ k {t-t') Pk {t')dt' (38) 

JO 

is equal in the two representations. Indeed, this arises from the fact that the memory 
function 7k(0> being the correlation function of the noise following FDT (p3[), must also 
be equal in the two representations. 

This will be only strictly valid if we could close exactly the hierarchy illustrated 
above, which clearly can also be seen as imposing constraints on the two sides of the 
density equations (Bl), or equivalently on the two representations of the noise. Since 
this is a far from practical program to carry on, we can try to identify those essential 
features of the system where consistency should be enforced, in order to render the whole 
approach simpler, but still satisfying the basic physical and mathematical properties. 

Consider now the two representations for the random noise ( |36|) and fl37|) . As noted 
above, the two should be equal to one another with respect to any averaging process we 
build. A limited expression of this constraint is expressed if we ensure the validity of 
the equalities for any order n, 

((flp-^t)) rtf®) = ((flp-^t)) vl n \t)} (39) 
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This means that the two sides of equation @ will be equal under projection onto density, 
density pairs, triplets and so on. It must be also noted that since these are all equal 
time averages they may be worked out exactly at any order and reduced to relations 
between the density correlation functions, and the potential. A clear example of this 
was given by the YBG equation (|i~9|), corresponding simply to the case n — 1 of the 
general condition (|3D|). 

A second example, that will have an important part in the following discussion, is 
represented by the case where n = 2, corresponding to the ensuring of the equality of 
the two random force representations under projection onto density pairs. Explicitly 
this can be written as, 

(P-pWP-cW^)) = (P- P (t)P-n(tW"\t)} (40) 

Of course, the contribution on the two sides of the integral term A k (t) (|38|) cancels 
out, because we are assuming in principle the exactness of the equality of the two 
representations at any level. Writing explicitly (f40|), we have, 



1 k 2 
~ k')wfc'(p_p(t)p-q(t)pk-k'(*)pk'(*)) - ^^c fc (p_p(t)p_ q (t)p k (t)) 

= (p_ p (t)p_ q (t)p k (t)) + n k (p_ p (t)p_ q (t) Pk (t)) (41) 
The terms (p_ p (£)p_ q (£)p k (£)) are easily calculated. Indeed, 

(p- p (t)p- q (t)p k (t)) = -[<p_ p (t)p_ q (t)p k (t)> + (p_ p (t)p_ q (t)p k (t)}] (42) 

and, 

p ■ k 

(p_ p (t)p_ q (t)p k (t)) = —-S^+^NSq 

(p_ p (t)p_ q (t)p k (t)) = V-±5 P+ ^NS P (43) 

Thus, rearranging equation (pTT|), considering also that (1 — ncujSk = 1, we finally have, 

1 k 2 
- TrZ!(k-k')^(p-p(t)p-(k-p)(t)Pk-k'(*)Pk'W) - -T(p_ P (t)p-(k-p)(*)Pk(*)> 
v k' P 

= --J}{ k2 - k • p nc p - k • (k - p) nc| k _ p |}S'pS'| k _p| (44) 

This represents indeed a further relation involving the bare potential and the static 
structure factor of the system, of the same type of the YBG equation (|T9|), but of higher 
order, now also involving the four density correlatione. We note that, if one consider 
the case p = in equation (|44"D, this simply reduces to YBG equation. 

To sum up, equations (^) for any n constitute a minimal set of exact equations 
for equilibrium correlation functions that must be satisfied if Newtons equations are to 
be satisfied. The first of such equations is known as the YBG equation, and has been 
widely studied. We have earlier shown that this consistency relationship has the effect 
of modelling the collective part of the particle motion in an optimal manner. The second 
consistency relationship has not yet been the subject of interest in liquid state theory. 
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It may be expected to play a role in ensuring that leading effects of phase relations 
between density variables are maintained, as well as some aspects of the dissipative 
effects of interactions that lead to decorrelation of the density variables. 

The set of equations ( |39| ) will ensure the fidelity of the Langevin process to the 
appropriate Newtonian dynamics on a single time slice in the history of the system. 
Indeed, they do not contain information about correlations between different time slices. 



6. Exact Calculation of the Memory Kernel 

Now we apply the FDT condition ( ^3|) to calculate the memory kernel. Following 
the definitions above, we have now two ways to proceed, by using the two equivalent 
representations for the random force rjk(t) (j36[) , ([37]) . 

Using the explicit form fl36|) , which contains in it the bare potential v k , it is 
straightforward to write an exact expression for the memory kernel applying FDT, 



ifW- — 



Nk 2 



k_ 

(3m 



2 \ 2 



/ 2 2 

In 



4 - l)NS k (t) + Q>. r,(0))V ik - r '(°)(k. r ro (t))V k - r »W) 

l,m 

+ -r^T,T, v k'Vk'ik-k')(-k-k , ')(p^ k ,io)p k/ io)p k ^(t) Pkl (t)) 
\ mV ) k' k" 

nk 2 { 

+ -^~ r c k \ 5> fc >(k ■ k')(p-k(0)p k -k'WPk'(*)> + 5> fc „(-k • k")(p-k-k»(0)Pk»(0)p k (t)> 

k' k" 

' k 2 • , „ 1 



(3m 2 V 

- f dt' lk (t-t') \ ^L n Nc k S k (t') + -L-J2v k <{-k- k')(p- k -k'(0)Pk'(0)Pk(t')) 
Jo am mv T7 



(3m 1 1 v v y 

Thus, to explicitly evaluate the memory kernel, it is necessary to introduce some kind 
of approximations to evaluate the multiple density correlations appearing in (|45|). 
On the contrary, when using the implicit representation (II) for the random force (|37|), 
the formally exact expression for the memory kernel, would be, 

li II] (t) = ^ {(p-k(O)pk(t)) + ^[(p-k(O)pk(t)) + (P-k(O)pk(t))] 

+ (fi fc ) 2 (p-k(0)p k (t)) +^V Tk (t-t')((p-k(0) + fijfeP-k(0))p k (t'))} (46) 

This contains, for example, correlations at unequal times of the second derivatives of 
densities which are not known, even in an approximate form. Thus, to pursue this 
second path, we must find a more suitable form for the random force to calculate its 
time correlation function. This is precisely the route that has been used in the original 
derivation of the mode coupling theory (MCT) |Tj|. Clearly this second possibility, 
despite the presence of these intractable objects, does not contain explicitly the bare 
potential. Of course, in an exact expression of the theory, we should have, 

7k /) W = 7 k //) W (47) 

but we caution that different types of approximations made on the two expressions for 
the memory kernel may lead to inconsistent results. 



(45) 
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7. Representation I: Random Phase Approximation and Mode Coupling 
Theory 

We now proceed to perform the simplest approximation to multiple density correlations, 
to solve the YBG equation ([19]) and at the same time calculate the memory function 
( fHf) in representation I. To do so, we consider the simplest case of a Gaussian random 
process ^k(^)- This implies, for example, that the fourth moment of the force can be 
decomposed in terms of second moments, as for example, 

(V-h(ti)V-h(t2)rju(t3)Vk(k)) = (V-k(t 1 )v--k(t 2 )}(Vk(h)Vk(t 4 )) 

+ (V-k(ti)rh(h))(V-k(h)r^(t A )) + (77-k(ti)?7k(t4)> fa-kfeO^feO) (48) 

Since we know that density variables are linear functions of the noise from (|3ll). 
also densities are then Gaussian variables. This means, for example, that the triplet 
correlation function appearing in the YBG equation ([19]) can be decomposed as, 

(P-k(*)Pk-k'(*)Pk'(*)) = (p-k(t))(pk-k>(t)p k ,(t)) + 

(Pk-kW>(P-k(£)Pk(t)> + (Pk'W>(P-k(i)Pk-k'(f)> (49) 
and by using the general result, 

(P-p(OPqW) = NS M [S{q) + NS^o] (50) 

we have, 

(p-k(t)pk-k^K' (t)> = a^m^ 

The result of inserting ( fyj) into (|i~9| ) is that only terms where k = k' give a 
finite contribution in the sum on the right hand side, but this term has been already 
isolated from the rest of the sum in (|i~9|) . Thus, this corresponds to a Random Phase 
Approximation (RPA) of the YBG equation. Therefore, we have that the structure 
factor is determined directly from the potential, via the relation, 

c k = ~Pv k (52) 

Now, consistent with this approximation, we can also calculate the density averages for 
unequal times in fl45|) , as for example, 

(p-k-k (O)pk' (0)pk-k" (t)f>k" (*)) = (P-k-k' (O)pk' (0))(Pk-v (t)Pk" (*)) 

+ (p-k-k'(0)pk-k''(*))(Pk'(0)pk»(*)> + (P-k-k' (0)pk»(*)>(Pk'(0)Pk-k»(*)> (53) 

Also, to complete the calculation of the memory kernel, there is a second minor 
approximation that we need to introduce to evaluate the second term in the right hand 
side of (|45|). Thus, we perform the average over the velocities, even if now we are not 
dealing with equal time averages, and so neglecting the correlations in the single particle 
kinetic energy, 

(EE( k • i-*(0))V^°>(k • r m (t))V k - r -«> » -f- 2 NS k (t) (54) 
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Thus, the Gaussian nature of the random process, and approximation fl5"3|), allow to 
rewrite (pEBj ) as, 

7k(t) = E K'(k • k') 2 + v k ,v k - k ,(k • k')(k • (k - k'))}S ]k ^\(t)S k ,(t) 

mv k k ,_^ k 

h^nP 1 ft 
+ -s—(c k + (3v k ) 2 S k (t) - n(c k + f3v k ) / dt'-Ykit - t')d t ,S k (t') (55) 
pm Jo 

At this point, we have to impose the YBG condition, consistent with Gaussian averages, 
that is the RPA condition d52|) , and this leads precisely to the well-known expression 
given by the Mode Coupling Theory (MCT) for the memory kernel, 

7k(t) A =V4iT7 E {( k • k ') 2 4' + ( k ■ k ')(k ■ (k - k'))c k ,c k _ k ,}5| k _ k1 (t)^(t). (56) 
prriK V k /^ k 

It is interesting to note that in the Gaussian approximation, by substituting the bare 
potential with the direct correlation function, the contributions to the memory kernel 
given by the integral term /q 7k(t — t')p k (t'), present in the noise are zero. If this were 
not the case, the expression for the memory kernel would not be explicit any more, but 
it would be instead a self-consistent relation. As it will be shown later, this will be 
also valid if one calculates consistently in a Gaussian approximation the memory kernel 
in representation II, projecting onto density pairs, but leaving aside the integral term. 



On the contrary, in the original derivation of MCT [0, using representation II, the 
approximations made are not all consistent with each other, and in that derivation, this 
term was a priori neglected, as explained in detail in section |9|. 

By introducing a factor 1/2 in (|56"D, and calling k' = p and k — k' = q, we can 
rewrite the memory kernel in the conventional MCT form, i.e. 

7xW MC = ^7 E + PCp) • e L (k)\%S q %(m q (t) (57) 

" q + p = k 
q^O 

where is the normalised density correlator, for which the equations of motion are 

now clearly, 

$k(t) + ^k$k(t) + /%k CT (* - t')^ k (t')dt' = (58) 
J o 

8. Projection onto Density Pairs 

As we have discussed earlier, the formal exact expression of the memory kernel 
in representation II (|46|) does not constitute a good starting point for making 
approximations and developing a simple theory, because of the presence of the second 
time derivatives of densities. A way to overcome this problem is to project the random 
force Q5TD onto the density pairs subspace, as it was done in the original derivation of 
MCT We, thus, leave unchanged the integral term A k (t) (PH[), which appears in the 
expression of the random force (|35"D, and in particular in (0), and project the remaining 
part, which is clearly equal to the minimum residual force !Fk(t) onto the density pairs 
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subspace. This choice follows naturally from the the YBG condition, which states that 
the residual force is orthogonal to the density at every time. 

For convenience, we define a reduced random force ffa(t) that is simply the random 
force from which we have subtracted the term A k (i) (P5|), such as, 

Vk(t) = rj k (t) - A k (t) (59) 

To be consistent, this reduced random force coincides with the residual force J-"k(t), but 
we prefer to use this notation at this stage to underline the fact that we are talking of 
a stochastic process, even if mapped onto a deterministic one. As for the pure random 



force, also nas of course two representations. Now, following |12| , we can rewrite 
the random force in representation II as, 

V*\t)= E (p-i(t)p- m (t)Ut) {II) )9(p,s,l,m)p p (t) Ps (t) + A k (t) (60) 

l,m,p,s 

The coefficients of the expansion (p-i(t)p- m (t)fj^(t)) can be easily calculated for the 
first two terms in the expression of the force (|37|). As we will show later, in the 
conventional MCT derivation, the integral term is simply neglected in the calculation 
of the memory kernel. The unknown matrix g(p, s, I, m) is the normalisation matrix of 
the projection. In these matrices the bare potential is now buried. These are defined 
implicitly via the relation, 

E(/'-p(^-sW/'iW/'mW)3(i,m,s'p') = * s ,^p ) p' (61) 

and it is clear from this that to explicitly evaluate g(p, s,l,m) there is the need to 
evaluate four densities correlation functions, and this can be done only by introducing 
some approximation for the stochastic process. For the sake of simplicity, in the following 
we neglect A k (£) but later we will show how it can easily be included in the derivation, 
and how it might be important to do so. 

The first important observation to make is that we can also think of performing the 
same projection onto density pairs as in ( |60D for the representation I of the random force 
(|36|). When this is done, it is easy to show that, at this level, this can be reproduced 
exactly, by making use of the definition of the normalisation matrix fl6T|), without the 
need to calculate it explicitly. Let us show this in detail. 

We want to calculate the coefficients of the projection for representation I of the 
reduced random force. We have, 

(P-l(t)P-xn(t)f]i (*)) = --^«C fe (p_i(t)p_ m (t)p k (t)) 

1 E( k - k ')^'(P-lWP-m(t)Pk-k'WPk'(t)) (62) 



mV k 



We note that these are static quantities, and thus the kinetic term in Q3q) can be 
exactly rewritten in terms of density when performing the averages, by separating out 
the velocity contribution. Thus, we insert these coefficients in the projection definition 
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([SO]), and we have, omitting the time dependence since everything is at the same time t, 
Vk ] = (P-iP-mViP)g(p, s, I, m)p p p s = -j^ nc k | ^2(p-iP- m pk)g{p, s, I, m) \ p p p s 



'.,m,p,s r p t s I l : m 

1 

k' P,s I l,m 



, llV 5Z( k ' k 'H' { H(/ ) -iP-mPk-k'Pk')fi , (p ) s > h m) \ ppp s (63) 



We can thus evaluate explicitly the sums enclosed by graphs in flBlf ) using (|6lD as, 

J2(P-iP-mPh)g{p, s, l,m)\ = ^p,o5 s ,k (64) 

J2(P-\P- m Pk-k>Pki)g(p, s, I, m) \ = 5 p ,k'^s,k-k' (65) 

l,m J 

and by inserting these results into (|6^) , remembering that pk=o = N, we find, 

k 2 1 

Vk\t) = -^nc k p k (t) - — Y.^ ■ k 'W'pk-k'(t)p k/ (t) (66) 

which is exactly the original expression we had for fj^\t) before projecting it. 

This result tells us that, since the two representations of the noise are intrinsically 
equal via the mapping between Newton and Langevin equations, also the projection 
of the reduced force in representation II onto density pairs in principle will be exact. 
Nonetheless, the problem with it lies in the fact that it cannot be calculated explicitly, 
because the definition (|6T| ) cannot be applied there to eliminate the unknown matrices 
g. Thus, differently from the original work on MCT [ Efl , where the projection onto 



pairs was just chosen as the simplest approximation and not investigated further, we 
can now say, through our knowledge of the Newtonian representation which lies beneath 
oyr stochastic process, that it is an exact step, if taken carefully, and it does not need 
further terms, like for example one could think to project onto density triplets and so 
on, as a better form of approximation. 

Nonetheless, this projection must be treated carefully. Indeed, it tells us that we 
have chosen the right subspace to project onto, but it also impose a further constraint 
on our system. As discussed previously, the equivalence between the two representations 
should, in principle, satisfy the constraints (|3^) at any order n. Here, since the use of 
representation II involves the projection onto density pairs, it is explicitly required that 
our system satisfies at least the second order type of constraint, expressed by fl4"4[). 

Thus, this equation must be taken into account when using representation II of the 
random force to calculate the memory kernel, because the projection onto density pairs 
is now only exact when this constraint is satisfied. 

8.1. Calculation of the Memory Kernel for the Projected Representation II 

We now proceed to calculate the memory function using the projection onto density 
pairs of representation II. Concerning the integral term Ak(t), we already know that 
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its contribution is null when we use the Gaussian approximation to evaluate the exact 
formula for the memory kernel in representation I In section 0, we will also 



reproduce the calculation of the kernel as in the original MCT derivation [12] and show 
that there the term A k (t) is simply neglected. Here, we will not neglect it and show 
that it is possible to incorporate it in the calculation, without rendering the calculation 
much harder. Thus, as illustrated in (BCT), we do not project it onto density pairs. 

Using the explicit expression for the coefficients of the projection for the second 
derivative term, given in (fy|) and fl43f) , we can write the 'exact' expression for the 
memory function, 'exact' because it is only to be considered exact when used in 
conjunction with the equation of constraint (@), as discussed in the previous paragraph. 
Thus, applying (|23|), and performing the sums that can be done explicitly so to isolate 
the various types of density correlators, we have, 

^ ® = ~Bmj^m, P ,J^ m ',p',s' [ fc2 - k ' m'nc m , - k ■ (k- m')nc fe _ m ,] [k 2 -k-mnc m - 
k-(k-m)nc k -m]Sm f Sk-rn'S m Sk-rn9(p, s,k-m,m)g(k-m',m',p', s')(p- p >(0)p-. s >(0)p p (t)p s (t)) 

- t. — ^rV] [k 2 -k-mnc m -k-(k-m)nc k - m }S m S k - m g(p,s, k - m,m)(p_ k (0)p p (t)p s (t)) 

pmbk rn,p,s 

- — — — V , , / [k 2 ~k-m'nc m/ -k-(k-m)nc k - m/ }S rn/ S k - ml g(k-m\m',p', s') 

(P-p'(0)p-s'(0)PkO)) - ^^ m , pls \k 2 ~k-mnc ml -k\k-m)nc k ^ ml }S 

^ k (*-O<P- P K%-v(0)p k (O> + -jAsSkit) + ^ fdt'^it - O(p-k(0)p k (t')) (67) 
o prnb^ NS k Jo 



This form has now the same structure as expression (|45|) for representation I. In other 
words, it is explicitly visible that it contains the time-dependent 4-point correlation 
(p_p/(0)p_ s '(0)pp(i)p s (£)), as well as triplets and pairs time-dependent correlations. 
Nonetheless, now the coupling coefficients do not contain the bare interactions, but 
already we can see that they are composed of static quantities, such as the direct 
correlation function, although clearly the unknown matrix g still needs to be treated. 

8.2. The Gaussian case 

In this section we want to apply the simple Gaussian approximation for the random 
process that we used to calculate the memory kernel in representation I, to the case of 



projected representation II. The first equation to study is of course fl44|) , because this is 
the one that guarantees the consistency of the whole scheme. Thus, recalling that for 
a gaussian noise multiple averages can be decomposed as in (|8|), the constraint (|4]) 
becomes in this simple case, 

1 k 2 

- —{k ■ pv p + k • (k - p)w| k -p|}iV 2 5'p5'|k-p| - -j{ S p,o + VkjA^S^ = 

- ^{k 2 - k • p nc p - k • (k - p) nc| k _ p |}S , p S'| k _p| (68) 
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Now using the condition c k = —f3v k , that we found applying the gaussian approximation 
to the YBG equation, we have cancellation of some terms, but then this would imply, 
{^p,o + $p,)i}NSk = S p S\k- p \ and this is clearly not an equality valid for all k and p, 
but only for the special cases k = p, or p = 0. Thus, we conclude that equation 
(fH|) cannot be satisfied by a simple gaussian approximation. This, indeed, implies 
that using representation II, projected onto density pairs, and proceeding with a pure 
gaussian approximation to calculate the memory kernel, the result found for it will 
not be the same as the one found with representation I (|56|), because the relation of 
constraint ( f44[) would be violated in this way. As a consequence of this, we have that, to 
use representation II, one must go a step further in the approximation of the noise than 
with representation I, where the simple gaussianity of the noise is sufficient to develop a 
consistent theory. This comes from the fact that representation I is intrinsically exact, 
while representation II is not, with respect to the projection onto pairs, and thus, using 
the latter, one must be more careful. 

To evaluate explicitly ([)?]), we need an expression for the normalisation matrix g. Using 
again gaussian properties of the noise to decompose density correlations, we have that 

(M) g ives > 

Using this result, and applying systematically gaussian decomposition of averages at 
equal and unequal times, we have the following result for the memory function, 

7k J) Gauss (t) = 2 pJ lk 2 N Y\ k2 - k • p ncp - k • (k - p) nc| k _ p |] 2 S p (t)S| k _ p |(t) 



k 



2 



nc k )S k (t) 



s 1 



7k(i-0^ fc (0 (70) 



PmS k 

The first thing to note here is that the contribution, given by the integral term, 
vanishes because [(1 — nc k ) — 1/Sk] =0, and this coincides with what we have found 
in representation I in gaussian approximation. Also, we can isolate from the sum the 
terms for which m is singular (m = and m = k), and, recalling equation (0), we can 
rewrite flTOf ) as, 

7k(t) = 7k(t) M + -5 * E^ 4 -2k-p^c p -2k-(k- P )n C | k _ p |]S p (t)5| k _ p |(t) (71) 

Thus, clearly, this way of proceeding does not reproduce MCT result and the reason of 
this lies in the fact that the two representations for the noise (|36|), ([37]) are not equivalent 
under projection onto density pairs using a simple Gaussian approximation, because, as 
shown above (|44D is not satisfied, except for the trivial wave-vectors cases. This means 
that further investigation of the theory must involve a new type of approximation for 
multiple density averages, that would satisfy firstly equation (f44l) . 



Dynamics of supercooled liquids: density fluctuations and Mode Coupling Theory 23 



8.3. The STLS Approximation: a step beyond Gaussianity 

We want to illustrate here a step in the direction of improving the gaussian 
approximation for density averages. Indeed, turning to consider other known 
approximations in literature, it may be possible to satisfy the constraint (144]) for a 
larger number of wave-vectors than the simple gaussian one. The example we will refer 
to is known as the Singwi-Tosi-Land-Sjolander approximation (STLS) |19[ . What we 
show here does not lead to a definitive result for corrections to MCT, but it is illustrative 
of the issues involved. 

We can consider now for example the triple average in the YBG equation (|19|). 
Recalling what happened with the Gaussian approximation, we can now think that in 
this new scheme, we will have, as dominant contributions to the sum, not only the 
Gaussian terms but also some corrections to them, that would be most meaningful if 
they in some way include the gaussian terms, so that the new scheme is actually an 
improvement of the old one. To say this mathematically, we can write explicitly, 

(p-k(t)Pk-k'(*)Pk'(*)> = N iKo + Ci}[S k > + N5v, ] 

+ N{8 kM , + C 2 }[S k + N5 K0 } + N{5 k ,, + C 3 }[S k + N5 KQ ] (72) 

where Ci,2,3 are the three corrections to the Gaussian terms, that are simply 5-f unctions. 
Thus, we are still assuming that the triplet correlation function decomposes into three 
separate contributions, as for the gaussian variables, but now these contributions will 
be more complicated. Following the literature |19| , the STLS consists in assuming the 
following, 

i(k-k')-r n (t) \p p ik'T„,(t) _ P *-r„,(t) V e *( k - k ')-(Mt)-r„,W) 



Pk-k'(*K'(*) = 5V (k_k } - r " (t) ^-y k '"v(*) = ^ e *-vW £ 



« ^ e *. V ( t )^ e i(k-k').(r n ( t )-r n ,(t)) ) = p k ( t )[S| k _ k ,| + NS^ V ] (73) 

ri n 

where the approximation made consists of replacing the sum over n on the second line 
of ( [T3"D by its average, as it was independent on n' . 

If we apply this to the triplet correlation function above, we have, 

(p-k(t)Pk-k'WPk'(t)) = N{5 K0 + S k }[S k , + N5 k , fi ] + 

+ iV{5 k , k ' + S\*- V \}[S k + N6 K0 ] + N{8 kl>0 + S k ,}[S k + N6 K0 ] (74) 

and the corrections £1,2,3 are simply the static structure factors, with wavenumber 
corresponding to the respective 5-function. A physical picture of what this represents 
is straightforward. Indeed, the gaussian decomposition of averages simply means that 
the densities at three points are only coupled two by two in all possible combinations, 
assuming that the third one is infinitely far from the others and thus, this is averaged out. 
The new STLS approach consists in accounting for some aspects of the correlations of 
the third particle. This is the reason why we now need two structure factors to describe 
our triplet correlations, while in the gaussian treatment we needed only one. 



Dynamics of supercooled liquids: density fluctuations and Mode Coupling Theory 24 



Now we apply the approximation (|7J]) to the YBG equation ([19]), thus obtaining a 
new solution for the direct correlation function of the system, i.e. 

c k = ~(3v k - JL £ «*(k • k')S(\k - k'|) (75) 

where again only the contribution driven by the k = k' 5-function survives, because 
the other two can be both shown to have spherical symmetry and thus, their integral is 
null. The result (|75|) clearly still contains the information of the RPA in the first term, 
but adds to it a new term to be taken into account for the calculation of S(k) from the 
potential. 

The next step now is to see what happens when applying the STLS approximation to 
the second equation of constraint (f4~4|). We have seen that the Gaussian approximation 
would satisfy this condition only for two singular wavevectors. Using the STLS 
approximation for the quadruplet and triplet correlations on the left hand side and 
([TED for the direct correlation function on the right hand side, (fHj) becomes 

/l v k'Qz- ' k')[SkS p Sk> + (S'ik'-pi + S'|k-k'-p|)»S'p»S'|k-p|] = 



k' 



(k- P )(k'-p) (k-(k-p))(k'-(k- P )) 

-£|k'-p| H 77j Ty2 >->|k-k'-p| 



5 p 5|k-p| (76) 



p* i— « (|k- P |) 2 

The first of the three terms on the left hand side gives no contribution due to spherical 
symmetry. Thus, to satisfy the constraint (f44|) in the STLS approximation, we have to 
satisfy the general condition for certain wave-vectors p, equivalently k — p, 

(k. k 'H < k - "><"'■■» (77) 

This implies the geometrical condition cos 9k,k> = cos 9k, p cos 9k', p , which is equivalent to 
sin 9 k tP sin 9 k', p cos 4>k,k' — 0. Being here k' arbitrary, the only general condition that can 
satisfy this constraint is that k and p are parallel, plus of course of both of them being 
zero. This is certainly an enlarged number of wave-vectors over which STLS satisfies 



the second order constraint ( 44]) with respect to the only two suitable values for the 
gaussian approximation. This type of scheme would possibly provide a better way of 
proceeding in calculation of the memory kernel, using the projection onto density pairs 
of representation II. A problem, in doing so, would be how to extend correctly the STLS 
approximation (1731) to unequal time correlations, and if then would then be possible to 
close the theory, by eliminating the bare potential, for example using the result of the 
YBG (0). 

These questions do not appear to have an easy answer, as well as the challenge 
of satisfying the constraint ([44]) for any wave-vector seems quite hard. That is why 
future work will be addressed to propose an alternative, in principle self-consistently 
determined, way to go beyond the gaussian closure, as well as to renormalize the bare 
interactions. 
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9. Comparison of our formalism with the original derivation of MCT 



Here, we want to link our formalism the traditional derivation of MCT [|12|], based on 
the projector operator formalism. In this way, we will be able to show in more detail 
the various steps that were originally taken, and compare them to our work. 

Thus, we analyze the original derivation, starting form when equations of motion 
are written for the density and current correlators, following the Zwanzig-Mori (ZM) 
formalism, where the variables which define the projection space are, indeed, the density 



Pk and the longitudinal current jk |12 |. In the subspace spanned by these variables, the 
projection operators are defined as V being the projector onto the subspace itself and 
Q = 1 — V the orthogonal projector. Thus, translating the notation in |TJj into our 
notation, we have the following equations of motion, 



V jk(f) J \n k /k o/VjkW/ Vo Dc(t-tf)J Vjk(t') J \R*(t) 

that combined together give equation (f2~T|), where R^(t) is the random force which was 
indicated as rj^t) in our notation. There are a few things to specify regarding (|78|) . 
Firstly, in the ZM formalism, the random force at the initial time t — 0, is defined for 
a general variable A as Ra(o) — QA, and so it is clear that Qp k = and the only 
non- vanishing random force is related to the current. Thus we have, 

flk(0) = Qjk(0) = -i/MPk(0) + fi fe pk(0)] (79) 

Also, in this formalism, we know that the random force evolves in time not with the 
Liouvillian operator C, but with a particular projection of it, and thus, 

R k (t) = exp (iQCQt)R k (0). (80) 

Nonetheless, in the MCT approach, the evolution of the random force is approximated 
with the Liouvillian evolution in time, i.e. with e tCt Rk(0). Thus, the approximation 
corresponds to neglect the difference between e tQCQt Rk(0) and e lCt R] i (0) i.e. 

Rk{t)*-^\pk(t) + tokPk(t)]. (81) 

In our derivation, using representation II, the true random force T)k(i) is given in (pJ7|). 
and it is clear that, apart from the factor — i/k that comes from using as a variable jk(^) 
instead of f>k(t), the two being related via the relation Pk(t) = ik]k(t), this is the same as 
the exact Rk(t) flSOl) . The comparison between ( |37|) and fl8lD allows us to give an exact 
expression for the difference between the true random force and the MCT approximated 
one. Indeed, it turns out to be the integral term Ak(i), defined in (|38D . We have shown 
how to include Ak(t) in the equations for the memory kernel. We have also shown that 
the contribution to the memory kernel coming from Ak(t) is strictly zero in RPA, in 
both representations. 

Thus, the starting point for the original derivation of the memory kernel in MCT, 
is equation (|67D, in which one neglects the integral term. The approximations that 
are then made are the following. Firstly, it is chosen to perform a kind of a gaussian 
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approximation for the four-point density average, which means that still the quadruple 
averages factorises into averages of pairs, but not all the possible combinations are 
considered, i.e. 

(p- P '(0)p-s'(0)Pp(£)Ps(t)> « (p- p <(0)p p (t))(p_ s ,(0) Ps (t))N 2 S p (t)S s (t)5 P:P ,5 s , s , (82) 

The same approach is followed to evaluate the normalization matrix g, using definition 
(|BT|), i.e. g(l,m,p,s) ~ o~L S o~m, P /N 2 S p S s . These types of simplified gaussian 
approximations can be substituted by the conventional ones that we have used, like 
for example (^), without changing the final result for the MCT memory kernel (|57|) , 
because of the symmetry that this has with respect to the wave-vectors p and q. 
Regarding, instead, the triple averages appearing in the coefficients, they can be exactly 
written as, for example, 

(p_i(t)p_ m (t)p k (t)) = N5i +mik SiS m S k (l + nc 3 ) (83) 

where C3 indicates the direct correlation function for triplets. Thus, based on the 
observation that for hard spheres systems it has been shown that the contribution 
given by C3 in (|83| ) is very small, this is neglected. Also, if one wanted to include 
them in the calculation, they can be easily calculated from simulations. Thus, this 
approximations for the triplet density correlations appears to be more advanced than the 
simple decomposition of it, generated by a gaussian noise choice, as in (f|9j). Nonetheless, 
this breaks the full consistency of the scheme of approximations in developing the theory. 

Eq. (^), but where the triplet direct correlation function is neglected, leads to a 
crucial simplification of the coefficients of the projection of the noise. Indeed, it implies 
a cancellation of the triplet term with the term in k 2 . Thus, one can see that this, in 
addition to the gaussian type of approximation for the quadruple correlations and the 
normalisation matrix discussed above, leads very easily to the well known expression for 
7^ CT (t) given in fl57|). 

Now let us analyze what the inclusion of the integral term A k (£) would produce 
in the present scheme. To calculate the full memory kernel we are left to evaluate the 
last term in (|6~7f) . Using the approximations just discussed, i.e. gaussian averages for 
4-point correlations and (|83|) for the triplets neglecting C3, we have that the last term 
reduces to, 

-T7T? ^2„ , , , i _ 7T-[^ 2 ~ k " m'nc m / - k • (k - m')nc k -m'}Sm'S q -m'5i' + m',k 

Nk 2 ^V,m',p',s' y f] m L 

+ ^^(P-k(0)pr(0)p m KO))U(/',m / ,p / ,/) ^W Tk (*-t / )(P-p'(0)p-s'(0)pk(t , )> 



(3m 



V[kWnc m ,+k-(k-m>c fc _ ro ,] [dt'^(t-t')(p_ m ,(0)p^ +m ,(0)p k (t')) (84) 
, Jo 



Thus, there is the need to calculate the dynamical average (p- m '(0)P-(km')(0)pk(^')) to 
evaluate the contribution of this term. We could use either the Gaussian approximation, 
as for the quadruple average (|8"2]), or a generalization for unequal times of the triplet 
expression fl83|). Nonetheless, in both cases, this contribution does not appear to be zero, 
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as in purely gaussian treatments, and moreover, it requires the introduction of another 
approximation, somewhat consistent with the others that have been made previously. 
In any case, it would be certainly different from zero, and its simple neglection, without 
a consistent justification of it, is clearly unsatisfying. 



Conclusions 



In this paper we have attempted to provide an alternative way to derive a mode-coupling 
theory. We have shown how Newton's equations can be represented by trial equations, 
and how these may be approximated. The approach starting from Newton's equations 
allows us to derive explicit and exact expressions for the the memory function and 
for the noise without the requirement of projection operators.lt thereby provides an 
alternative view that may lead to new insights. The resulting equations, solved under 
the assumption that the noise, and consequently the density fluctuations of the liquid, 
is gaussian distributed, are equivalent to the random-phase-approximation for the static 
structure factor and to the well known ideal mode coupling theory (MCT) equations 
for the dynamics. This suggests that MCT is a theory of fluid dynamics that becomes 
exact in a mean- field limit. This possibility was suggested some time ago on the basis 
of the analogies between the equations describing the schematic MCT models and the 
the dynamics of the order parameter in disordered p-spin models, solved under strict 



mean-field approximation |TJL IT 



We have also investigated in detail the basic meaning of the less controlled steps 



intrinsic in the derivation of conventional MCT [TJ3]. One important point is the 
demonstration that, if the projection of the noise over the density pairs could be carried 
out without approximations, the resulting expression would be exact. The evaluation 
of the quadruplet correlation functions under gaussian statistics, an approximation 
employed in the evaluation of the memory function when projected in the density 
pair subspace constitutes, in conventional MCT, the step equivalent to the RPA 
approximation. Hence, our approach suggests that the conventional MCT is exact in 
the RPA limit. This offers the possibility of improving the theory by improving the 
approximation of the four-point correlation functions in the same way as one improves 
the static structure factors on moving beyond RPA. We have illustrated an example of 
this by using the STLS approximation (|73|), which leads to a different closure for the 
YBG equation (|75|), certainly more advanced than the simple RPA. 

The implication of this is clearly not that ideal MCT is a theory which can be 
applied only in the weak-coupling limit, where RPA becomes exact. In fact, the RPA 
limit that we have studied may in principle be corrected in two generic manners. Firstly, 
the treatment of the static correlations may be improved by use of a better direct 
correlation function than that implied by the consistent theory. On the other hand, 
Kawasaki has presented an interesting derivation of MCT based on the (quadratic) 
density functional Ramakrishnan-Yussouf free energy of a liquid, where the effective 
interaction between the density pairs is exactly —c(k)/(3, thereby accomplishing the 
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limited form of renormalisation implied by RPA. In addressing this issue we note that 
the comparison between the MCT prediction and hard-core type experimental ^ 
or numerical [28| results in dense liquid states are sometimes astonishingly good, even 
for network forming liquids |29 , |30| . 

However, the presence of the bare interactions in our formulation of the theory 
constitutes a strong limitation to our approach, and, to our present knowledge, the 
gaussian approximation, both for static and dynamics, appears to be the only immediate 
fully consistent closure. Thus, the main problem in developing further this formalism to 
better types of approximations is how to renormalise more effectively the interactions, 
othen than simply RPA. We will attempt to do so in future work. 
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